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Abstract 

The following inverse problem is discussed. A static electromagnetic field generated by a limited system 
of charges and currents is supposed to be known with its first derivatives at a point somewhere far from the 
system. This allows to reconstruct the position of the system, its net charge, and the electric and magnetic 
moments of the system. 



1 Statement of the problem. 

By the inverse problem of electrodynamics we mean the problem of reconstructing the charge and current 
densities from the known electromagnetic field they create. Obviously, the solution of the inverse problem of 
electrodynamics depends substantially on the presence of materials, the region in which the electromagnetic 
field is specified, ets. It is well known that the value of an analytic function and all its derivatives at some point 
allow to reconstruct the value of the function in the domain of definition by use of the Taylor series. Applied 
to the electrostatic inverse problem it means that if we know the value of electromagnetic field in some point 
and all its derivatives at this point, we can reconstruct the field in the whole space. Furthermore, the Maxwell 
equations allow to calculate the charge and current density functions p(r) and j{r) respectively: 
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Pir) ^ ^divEir) j(r) = —rot i/(r). 

Here E{r) is the electric field and H{r) is the magnetic field. 

In actual practice we are able to measure the field and only few first derivatives with an inevitable error. This 
means that we can reconstruct the field only in a small vicinity around the point where the field is measured. 
But the real problem encountered in practice is to calculate the field and charge distribution far from the point 
where the field is measured. In order to do it we have to know the value of field and its derivatives up to a very 
high order and with high accuracy. We present here another approach based on the multipole expansion of the 
field far from the static collection of charges and currents. 

Suppose we know the value of electric E and magnetic H fields at some point far from limited system of 
charges and currents. We know also the first derivatives of fields Eij = dEi/dxj and Hij = dHi/dxj. There is 
to define the charge q, the electric d and magnetic m dipole moments of the distant electromagnetic system, 
and the location r of the system. Let the coordinate origin be at the point where the fields are defined. We 
start with the well known formula for the fields of electric and magnetic dipoles 00 
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Taking the derivatives of these expressions with respect to coordinates we find 
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where Sij is the Kronecker symbol. 


One can see from Eqs. 


111.31 and 11.411 that the tensors Eij 


and Hij are 



symmetricaL This follows also from the Maxwell equations, namely from rot £? = and rot H — 0. Besides, 
the Maxwell equations div E = and div H = give 

Ell + E22 + Ess ^ Hii + H22 + Hm = (1.5) 

We consider expressions II. in - 11.41 as a set of equations in 10 scalar unknowns q, r, d and m. But this set 
consist of 16 equations with regard to the conditions 11.51 . and symmetry of the tensors Eij and Hij. This 
allows us to formulate the inverse problem in more general sence. Namely, we can consider fields E and H as 
being originated by difi^erent sources. Let us denote by the radius-vector of a charge collection with a total 
charge q and an electric dipole moment d, and by rm the radius- vector of a system of currents with a magnetic 
dipole moment m. Then the set of equations II. in - 11.411 splits into two independent sets for and rm- In 
solving these equations we can obtain in particular r, = r™,, which means that the electric field E and magnetic 
field H are generated by the same electromagnetic system. 

It should be pointed out, that in calculation of fields far from an electromagnetic system we can neglect the 
higher terms of multipole expansion. The accuracy of such representation depends on the ratio between the 
sizes of the system and the distance between the system and observer. In solving the inverse problem we find 
only the distance |r|, but not the sizes of system. In order to estimate the accuracy of the received solution one 
have to calculate the sizes of the system by some independent method. For example one can solve the inverse 
problem for a few different positions of the observer and then one can estimate the sizes of the system. One 
can also investigate the field in the vicinity of the observer in order to find out whether this field is of dipole 
nature. These problems are not considered in this paper. In order to exclude such questions we consider further 
a point-like charge with an electric dipole moment and a point-like magnetic moment. In the next two sections 
we solve the inverse problem first for the equations 11. H and 11.31 and then for equations 11.21 and 11.41 . 



2 Inverse problem for a charge and electric dipole. 

Let us solve the equations 11.11 and 11.31 with respect to the unknown g, d and r^ . We choose the coordinate 
system as follows: the coordinate origin is placed at the point where the field is specified, the x axis is directed 
along the vector E and the y axis is aligned with the principal normal to the electric field line. The vector of 
the principal normal n is defined by the equality 

- 

^ kE~ds 

where k is the curvature of the field line and ds is the displacement along the field line. Taking derivative from 
E we find 

n^^[E[DE]] 

where D — (iSV)iS with the projections Di — EkEik (a summation over repeated indices is implied). The 
projections of the vector n on the axes of an arbitrary orthogonal coordinate system are expressed in terms of 

rii = -£^EkEj{EjEik — EiEkj). 

If the X axis is directed in the sense of the vector E, then 

n = -^{0, E12, Eis}. (2.1) 

With the y axis directed in the sense of the vector n we have E13 = and £'12 > 0. It can be seen from Eq. 
11.11 that the vectors E, d and r lie in the same plane; therefore, the problem is reduced to a two-dimensional 



^We omit the subscripts q and m having in mind that in the section 2 we are dealing only with rq and in section 3 only with r„ 
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one. In the x, y, z coordinate system we have E — {E,0,0), d — {dx,dy,0), and ,£23 = as a consequence of 
the axial symmetry of the field. Thus, five functions of the coordinates E^, Ey, En, E12 and E33 are known. 
They can be used to find five unknowns q, d^, dy, ri and r2. Egs. ll.H and 11.31 in the {x,y) plane take the 
form 
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The charge q and the components of the dipole moment can be easily expressed from Eqs. (12.211 . 12.311 and 12.61 : 
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dy 
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q = -2r'^£'33 - SxrE 
Substituting these formulas into Eqs. 12.411 and 12. 5L we find 
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Taking into account that r — ^a;^ + y^, we get two equations depending only on x and y. Thus, the problem 
is reduced to the solution of the last system of equations 12.1UI I. By simple algebraic manipulations this system 
can be reduced to the following form (r 7^ 0): 



'(2£;33 + Sii) + 2/'(£;3 



Eii) + 2yxEi2 = 
0. 



-xy{E33 + 2Eii) - y^Ei2 + x^Ei2 + ZyE 
Let us introduce the designations 

El — E22 — E33, E2 = i533 — -Ell, F3 — Ell — E22, 
S = ' 
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Ef2 + FiF2. 

Then the solution of Eq. H2.11^ can be written in the form 
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In order to calculate the dipole moment, we find the squared radius- vector r. 
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S^{4E!, + Ei: 
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Substituting Eqs. I)2.12|l . I|2.13|l and 12.14^ into Eq. 1)2. 7|l . we find the dipole moment 

= tE'°{Ei2{E? + {-Ei2±Sf)T2E?S) 
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"(5)4(F2 + (-£12 ±5)2)5/2 

The sign ± in the formulas indicates the existence of two solutions of the initial system of equations 12.212.61 . 
Physically this means that the same field with its derivatives may be created at the given point by two different 
sources located at different places. One can find some specific cases of this solution in Ref. 
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3 Inverse problem for a magnetic dipole moment. 



In order to find the position vector r and tlie magnetic moment m of a particle generating the field H we solve 
equations 11.211 and 11.41 . Now we align the axis x along the vector H and the axis y along the principal normal 
to the magnetic field line. Repeating the reasoning of the previous section we get 

= ^ir^)y^-^'^y (3.2) 

H,, = „3(^-6^ + 15.^(^ (3.3) 

H,, = -3 "^-" + '"^" +15:.yi^ (3.4) 

H,, = -si^ (3.5) 

Thus, we have 5 equations for 4 unknowns. Hence, the system of equations (13.111 - 13.51 is a overdetermined 
one. It means that the components of H and Hij are not independent. The fact that the field H is produced 
by a magnetic moment places a constraint on H and Hij . The corresponding equation will be found latter (see 
Eq. lISIIHJ)- Substituting (rm) from Eq. ifT^t into Eqs. (glU - IpOt . we find 

H = ^H33X-^ (3.6) 

= -ff33y-^ (3.7) 

u u a^xX 2H33 

Hu = H33 - 6 — — - (3.8) 

TT m^y + myX H33 

H12 = ~3 bxy—- (3.9) 



Eliminating rrix and my we get 



= Haa+x^^ ~ex^ (3.10) 
H12 = sjy + xy:^ (3.11) 



and for the square of the distance 
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We suppose here that Hzd. 7^ 0. The case = will be considered latter. Using Eqs. I3.1U1 - 13.111 one can 
find 

3^(3^33 + H.) (3^3^ 
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Now we express r from Eq. iTTIi 
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Let us consider the specific case H^^ — 0. It follows from Eq. 13.511 that (rm) — 0, which gives immediately 
a; = 0. Eqs. (13.111 - 13.41 give the whole solution in this case 



y = 3— = -27—3- my = (3.17) 



Thus, Eqs. H3.13|l - H3.16|l and in particular case Eq. I)3.17|l give the solution of the problem. 

It was mentioned above that the system of Eqs. 13.111 - 13.511 is a overdetermined one. Let us find relation 
between H and Hij. On the one hand it follows from Eq. I)3.12|l that 

2 9H^{3Hs3 + Hii) . , 
^ 4^1 ' 

On the other hand relations 13.131 and 13.141 give 



This gives the desired equation 



In case -ff33 = we have from Eqs. (gllj - ^iB Hii = 0. 

If the condition 13.181 is not fulfilled, the field H is not produced by a magnetic dipole. 
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